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Abstract 

^A  new  method  is  presented  for  determining  the  moments  of  nuclear 
magnetic  resonance  absorption  lines  from  the  shape  of  either  the  free 
Induction  decay  or  that  of  the  echo.  Unlike  previously  used  techniques, 
this  method  does  not  require  the  assumption  of  an  analytical  function  for 
the  lineshape  or  the  fitting  of  the  experimental  decay  with  a  polynomial. 

A  fast,  suitably  precise  and  numerically  stable  algorithm  has  been  developed 
for  performing  the  integration  required  by  the  new  method. 
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The  moments  of  the  NMR  absorption  line  of  a  solid  are  important 
structure-sensitive  parameters.  The  increasing  application  of  NMR  rf- 
pulse  methods  to  solids  has  led  to  an  interest  in  determining  the  moments, 
especially  the  second  moment,  from  the  data  obtained  in  such  experiments. 

Two  different  approaches  have  been  employed  to  extract  the  moments  from 
the  shape  of  either  the  free  induction  decay  (FID)  or  that  of  the  solid 
echo. 

In  one,  the  experimental  decay  curve  is  fitted  by  an  assumed  analytical 
function1  and  the  moments  are  then  calculated  from  the  parameters  describ¬ 
ing  the  optimal  fit.  The  main  drawback  of  this  method  is  the  required 
assumption  of  a  particular  absorption  lineshape. 

In  the  second  approach,  the  beginning  of  the  experimental  decay  curve 

2n 

Is  fitted  by  a  polynomial.2  The  coefficient  of  t  in  this  expansion  Is 
then  set  equal  to  (-1)  M2k/(2k.'),  where  Mn  denotes  the  nth  moment.  The 
odd  moments  are  obtainable  in  the  same  way  from  the  out-of-phase  component 
of  the  induced  signal.  While  avoiding  any  assumption  about  the  lineshape, 
this  method  is  unreliable  from  a  numerical  point  of  view  since  very  good 
fits  may  often  be  obtained  with  very  different  polynomials. 

In  this  note  we  present  a  new  method  for  determining  the  moments  which 
requires  neither  an  assumption  about  the  functional  form  of  the  lineshape 
nor  the  fitting  of  the  experimental  decay  by  an  analytical  funct ion. 

Theoretical 

The  nth  moment  of  the  absorption  line  g(tw)  is  defined  as 
4®  "to 

Mn  =  J  w  g(u)0-Hi))duu  /  J  g(u>)du>  , 

-00  -00 

where  uj  is  the  center  of  the  band,  chosen  such  that 


(1) 
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M,  =  0  .  (2) 

The  shape  of  the  FID  (or  of  the  echo)  obtained  by  phase  sensitive  detec¬ 
tion  is 

■tto 

G(t,a)  =  Re  {exp[  1  ( Qt-fot) ]  J  g(uu)exp(  - iuut) duo}  ,  (3) 

-CD 

where  Q  is  the  irradiation  frequency  and  a  is  the  phase  of  the  detector. 

If  a  pair  of  mutually  orthogonal  phase  sensitive  detectors  is  used,  for 
the  first  of  which  a  =  cp,  one  obtains  two  FID  curves  described  by 

4oo 

G  (t)  =  G(t,cp)  =  Re{exp[i(nt-Kp)]  J  g(uj)exp(-iuut)duu}  , 

^  -00 

j 

G cp'(t)  =G(t,cp-f~u)  =  -Im(exp[  i  ( Ot-Kp)  ]  J  g(uu)exp( -itut)duu}  .  (4) 

—00 

The  functions  G^(t)  and  G^'ft)  satisfy  the  relations 
G^-t)  =  GCp(t)  C°S2CP  -  G^'Ct)  sin2cp 

V(-*)  =  sin2cp  -  G^'Ct)  cos2cp  .  (5) 

For  t=0  these  relations  lead  to  the  condition 

tarfP  =  -G^OJ/G^O)  (6) 

which  provides  an  operational  definition  of  the  phase  cp.  In  particular, 

It  makes  it  possible  to  adjust  the  phase  sensitive  detector5  until  G^'(0)=0 
in  which  case  <p=0  ( if  G^O)  >0). 

Assuming  that  the  phase  has  been  properly  adjusted,  we  introduce  the 
notation 

G||(t)  =  G(t,0)  and  G|(  t)  =  G(  t.Tt/2)  .  (7) 


Equation  (5)  now  reduces  to 


and 
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G |  j ( -t)  =  G| | ( t) 


Gjj-t)  =  -Gjj  t) 


(8) 


From  Eq.  (4)  it  follows  that 

+03 

g(u)Q+<u)  =  (  1/2tt)  J  [G||(t)  -  iGj^(t)]exp[i(w-6)t]dt  ,  (9) 

where  6  =  fl  -  ro  . 

o 

By  substituting  Eq.  (9)  in  the  definition  of  the  first  moment  and 
using  Eq.  (8)  to  simplify  the  result,  we  obtain 


M|  =6  +  Gj^(  0)/G|  |  (  0)  ,  (10) 

where  the  dot  denotes  the  derivative  with  respect  to  .time.  By  Eq.  (2), 
this  implies  that 

6  =  fl  -0)o  =  -Gj_(0)/G,|(0)  (11) 

which  may  serve  as  an  operational  definition  of  6.  Accordingly,  to 
satisfy  experimentally  the  on-resonance  condition  6=0,  one  simply  adjusts 
U)q  so  that  G^(0)=0,  with  G|j(0)^0  . 

With  this  adjustment  made,  we  may  shift  the  origin  of  the  frequency 
scale  by  putting  U)Q=0.  Equation  (9)  then  becomes 

g(oj)  =(l/2rr)  J*  [G|  |  r(  t)  -  lG^r(t)]exp(  itot) dt  ,  (12) 

where  the  superscript  r  refers  to  the  on-resonance  condition.  The 
even  moments  of  g(ty)  are  identical  with  the  corresponding  moments  of  its 
symmetrical  part  which,  by  Eqs.  (12)  and  (8),  can  be  written  as 

gs(u?)  =  (  1/tt)  Re  J*  G.  |  r(  t)  exp(  icut)dt  .  (13) 

o  *  ' 

Correspondingly,  the  odd  moments  of  g(tu)  are  the  moments  of  its  anti- 
symmetrical  part 


CO 

9  (ti>)  =  (  1/tt)  Im  J  Gi  r(  t)exp(  icut)dt  . 

0  -*■ 

Note  that  Gj^r(t)=0  if  the  absorption  line  is  symmetrical. 
Next,  we  introduce  quantities  Qn(  e)  defined  as 


-ko  n 

(*  U) 


*"<e)  =  i-^^  9Wd”  - 


where  e  Is  a  positive  parameter  which  has  the  dimension  of  time,  and  v 
equals  1  for  even  n  and  0  for  odd  n.  The  quantity  Qn( e)  is  connected  to 
the  nth  moment  by  the  extrapolation  formula 


Mn  =  Mm  <in(e)/G,,r(0)  . 

By  using  Eqs.  (13) -(15)  one  may  write 


<We)  =  1  G||^^f2n(e*t)dt 


(17a) 


^n+l^  I  GJ_  ^f2n+l^e,t^dt 


(17b) 


in  which  the  functions  f n( e, t )  are  given  as 


f2n<  =  (,/n)J 


exp(  iu)t )  duu 


f2n+l(6,,)  =  '(i/n,i  t-^cuu)2"*2 


exp(  iuut )  dto 


(18a) 


(18b) 


For  e  A  0  these  integrals  converge  and  their  values  can  be  calculated 
directly  from  the  residue  theorem.  The  resulting  formulas  for  the  first 
several  functions  f n( e, t )  are  (  for  positive  t)  : 


<t<rr- ;  »n***iw> 
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f  i(  e,  t)  =  ( 1/ e2)exp(  -  |)  , 

fzU.t)  =  ( 1/e3)exp(  -  ~)  cos(^|r  +  J)  , 

fs(  e, t)  =  ( l/e4)exp(  -  ~)  , 

f*U,t)  =  ( 1/3e5){exp(  -  “)  +  2  exp(  -  cos(^-  +  |)}  , 

f s(  e, t)  =  ( l/3e6){exp(-  |)  +  2  exp(  -  ~)  cos(^||)}  .  (19) 

The  individual  moments  can  now  be  obtained  from  the  observed  decay 
functions  G|  j  r(  t)  and  Gjr(t)  by  combining  Eqs.  (16),  (17),  and  (19).  It 
Is  important  to  note  that  the  limit,  Eq.  (16),  is  reached  in  each  case 
with  a  zero  slope  since,  from  Eq.  (15), 

lim  dQ,  (e)/de  =  0  .  (20) 

e-»0  n  ' 

This  greatly  facilitates  the  extrapolation  of  Qn(  e)  for  s-»0. 

The  off-resonance  condition  of  Eq.  (11)  may  lead  to  experimental 
error  when  either  the  dead  time  or  the  experimental  noise  limits  the  data 
available  in  the  vicinity  of  t=0.  For  this  reason,  we  now  consider  the 
effects  upon  the  moments  of  having  6^0.  It  may  be  shown  that 

,|^-Q 

6=  lim  J  Gj_(  t)  f ,(  e,  t)  dt/Gj  j  ( 0)  .  (21) 

Therefore,  the  extrapolation  procedure  used  to  determine  the  moments  can 
serve  also  for  an  accurate  determination  of  6.  By  comparing  Eq.  (9)  with 
Eq.  (12)  one  obtains  the  relations 

G|  |  P(  t)  =  Gj  j  ( t)  cos(6t)  -  G^(t)  s!n(6t) 

Gj^r(t)  =Gj(t)  cos(  6t)  +  Gj  j  ( t)  sin(6t)  (22) 

which,  once  6  is  known,  enables  one  to  calculate  the  moments  from  the  off- 


resonance  decays  Gii(t)  and  Gi(t). 


Off-resonance  measurements  have  lower 
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sensitivity  than  on-resonance  observations  and  there  are  additional  ex¬ 
perimental  errors  due  to  the  oscillatory  behavior  of  G|  |  ( t)  and  G^(t). 
Therefore,  Eqs.  (21)  and  (22)  are  best  used  to  correct  for  small  devia¬ 
tions  from  resonance. 

Comments  on  the  Numerical  Procedure 

For  n  >  2  and  for  small  e,  the  functions  fn(e,t)  rapidly  oscillate. 
This  leads  to  substantial  problems  in  the  numerical  integration  required 
by  Eq.  (17).  The  simple  Simpson  rule  proved,  in  fact,  to  be  inadequate 
with  any  experimentally  reasonable  spacing  between  the  digitized  data. 

The  choice  of  a  smaller  integration  step,  combined  with  suitable  interpola¬ 
tion  of  the  experimental  data,  increases  the  numerical  stability  but  the 
calculation  time  increases  as  well.  An  algorithm  was  developed  which 
proved  to  be  both  fast  and  numerically  stable. 

A  quadratic  Interpolation  polynomial  y j ( t)  =  ajt2  +  bjt  +  Cj  is 
determined  for  each  set  of  three  consecutive  experimental  points  G(t._j), 
G(tj),  and  G(tj+j).*  The  contribution  to  Qn(  e)  arising  from  the  interval 
Tl  ■<2<tl-l+tl>'2<tl+‘l+l>>  is  calculated  as 

qn!(e)  =  J*  G(t)fn(e,t)dt  =  y5( t) fp( e, t)dt  ,  (23) 

using  explicit  formulas  for  the  integral  on  the  right-hand  side  of  this 
equation.  The  quantity  Qn(  e)  is  then  set  equal  to  the  sum  of  all  qn'(e) 
plus  a  similarly  obtained  correction  term  arising  from  the  initial  interval 
<0,|( t i+t2)>,  where  t|  is  the  coordinate  of  the  first  experimental  point. 

At  this  point  we  have  to  stress  that  the  difficulties  connected  both 
with  the  numerical  stability  and  with  the  sensitivity  of  the  results  to 
experimental  errors  increase  drastically  with  increasing  order  of  the 
moment.  We  feel  that  with  reasonably  good  experimental  data  it  is  possible 
to  obtain  at  most  reliable  second  and  third  moments  and  semiquant i tat ive 
information  about  the  fourth  moment. 
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Results  and  Discussion 

In  order  to  test  the  proposed  method,  we  have  performed  the  numerical 
integration  required  in  Eq.  (17a)  for  H2,  using  for  G| | r( t)  analytically 
generated  sets  of  "data"  as  well  as  experimentally  obtained  digitized 
data.  Gaussian  and  more  complicated  decay  curves1  were  generated  with 
durations  corresponding  to  second  moments  ranging  from  1  to  16  Gauss2 
for  fluorine.  The  G | |  r( t )  values  were  calculated  with  a  spacing  of  1  nsec, 
which  is  experimentally  feasible  with  a  fast  digitizer.  In  Fig.  1  we  give 
the  results  obtained  for  the  ratio  0z(  e)/M2  as  a  function  of  e  for  several 
of  the  generated  decays. 

The  values  obtained  for  the  second  and  higher  moments  are  very  sensi¬ 
tive  to  the  shape  of  the  induction  signal4  G( t)  at  small  values  of  t.  The 
quantity  Q.n(  e)  depends  upon  the  full  range  of  G(t)  but  the  weight  given 
to  G(t)  for  large  values  of  | t|  decreases  rapidly  with  decreasing  e.  This 
does  not  present  any  problem  If  the  shape  of  G(t)  is  known  for  all  values, 
of  time,  as  in  the  case  of  echo  experiments  on  sol  ids— either  the  regular 
echo2  or  the  recently  reported  "magic"  echo. 5  However,  the  situation  Is 
substantially  different  for  the  free  induction  decay,  for  which  data  are 
not  obtainable  during  the  dead  time  cO.t^  following  the  rf  pulse. 

The  latter  point  is  shown  in  Fig.  2  which  reproduces  the  experimental, 
time  averaged  FID  of  the  fluorine  NMR  in  solid  KAsF*  powder  at  room  tem¬ 
perature,  observed  with  a  pulse  spectrometer  operating  at  25  MHz.  The 
data  are  digitized  for  every  microsecond  and  the  dead  time  is  about  12 
p.sec.  This  dead  time  leads  to  the  drastic  divergence  in  the  values  obtained 
for  Q2(  s)  at  small  values  of  e,  illustrated  in  Fig.  1  by  the  numerical 
results  corresponding  to  the  experimental  FID  in  Fig.  2.  If  the  divergence 
occurs  before  0.2(e)  has  approached  M2  to  a  suitable  degree  of  accuracy, 
the  data  are  inadequate  for  the  determination.  It  may  be  shown  that  the 


Values  of  Q2( e)  calculated  numerically  as  a  function  of  e 
for  several  induction  decay  curves  corresponding  to  different 
values  of  M2  for  fluorine  nuclei.  The  solid  lines  are  for 

"data"  points  generated  with  G| | r( t)  =  exp( -aztz/2) ;  the  dashed 

line,  withG||r(t)  =  exp(  -aztz/2)  x  ( l/bt)sin  bt.  The  points 

with  dashed  line  are  for  the  experimentally  obtained  G  t  j  r(  t ) 
given  in  Fig.  2.  II 
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Fig.  2.  The  ttme  averaged  l9F  free  induction  decay  observed  at  25  MHz 
for  KAsF6  powder  at  room  temperature.  The  data  are  digitized 
at  I  psec  intervals  following  a  dead  time  of  about  12  psec. 

Only  the  In-phase  component  Giir(t)  of  the  induction  signal 
was  observed. 
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fault  lies  in  the  lack  of  data  and  not  in  the  mathematical  procedure 
employed.  In  a  case  such  as  this,  the  determination  of  M2  by  fitting 
the  decay  data  with  a  series  expansion2  would  give  what  appears  to  be 
a  mathematically  definite  answer  but  which  is  physically  unreliable 

On  physical  grounds,  one  would  expect  the  divergence  to  occur  when 
c  <  tj,  which  is  verified  by  our  numerical  results,  those  in  Fig.  1  for 
KAsF6  being  typical.  It  is  seen  that  Q2(  e)  diverges  for  e  <  6  Usee, 
with  td  about  12  g-sec.  Moreover,  a  smaller  value  of  e  is  required  for 
Q,2(  g)  to  approach  a  larger  value  of  M2.  This  enables  limits  to  be  placed 
upon  td  for  establishing  M2  to  a  given  accuracy.  Our  results  in  Fig.  1 
indicate  that  if  an  accuracy  of  better  than  5  percent  is  desired  in  M2, 
the  dead  time  must  meet  the  condition 

td  <  2. 5x10  V  /mT  ,  (24) 

where  t  .  is  in  Usee  and  M2  Is  in  rad2sec”2. 
a  • 

Alternately,  one  can  attempt  to  correct  the  beginning  of  the  decay 
curve  for  the  distortion  due  to  slow  recovery  of  the  receiver  from  satura¬ 
tion.  A  method  has  been  described  by  Barnaal  and  Lowe6  for  this  purpose. 
Of  course  such  recovery  transients  are  avoided  in  the  echo  technique,  which 
is  an  advantage  of  the  echo  as  compared  to  the  FID. 
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